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ABSTRACT 


A  preliminary  theoretical  investigation  of  the  propagation 
of  a  quasi-static  E  mode  in  a  plasma  column  with  non-uniform  electron 
density  distribution  in  the  cross-section,  has  already  been  performed 
in  the  special  case  of  a  density  independent  £5/  ratio  (Technical 

Summary  Report  R.2). 

Starting  from  Maxwell's  equations,  a  more  general  case  is 
here  considered,  where  a  parabolic  distribution  has  been_,taken  as 
an  analytical  model  for  the  electron  density  profile  in  the  cross  - 
section. 

Electron  density  is  supposed  to  vanish  at  a  point  outside 
the  plasma  column  ;  in  the  limit  for  the  zero  point  to  go  to  infinity, 
uniform  case  equations  are  obtained. 

Letting  the  metal  walls  move  to  infinity,  or  coincide  with 
the  plasma  boundary,  the  completely  filled  waveguide  and  the  plasma 
column  in  free  space  dispersion  curves  have  been  verified. 

The  Ritz  variational  method  is  used  to  minimize  the  error 
over  the  cross  section  when  an  approximation  function  for  the  field 


distribution  in  the  plasma  is  assumed. 
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1 .  Introdution 

For  analitical  simplicity  let  us  consider  a  plane  geometry 
instead  of  a  cylindrical  one  and  let  us  suppose  that  a  slab  of 
inhomogeneous  gyrotropic  plasma  is  located  between  two  parallel  lossless 
metal  plates,  the  distance  between  them  being  2D  and  the  slab  thickness 
being  2d. 

The  coordinate  reference  frame  is  chosen  so  that  the  slab 
is  perpendicular  to  the  x  axis  with  its  surfaces  at  +  d,  and  the  z 
axis  in  the  direction  of  propagating  waves. 

The  magnetic  field  B  points  along  the  plates  in  the  direction 
of  the  z  axis.  The  plasma  slab  is  inhomogeneous  along  the  x  axis  and 
infinite  in  the  indirection. 

The  field  components  and  the  electron  density  distribution 
function  are  expressed  in  terms  of  the  independent  variable 
being  the  longitudinal  propagation  constant. 

The  whole  structure  and  the  cross  sectional  distributions 
assumed  for  the  E^  component  and  for  the  electron  density  are  shown 
in  f ig.  1 . 

The  analysis  proceeds  by  calculating  the  longitudinal 
components  of  the  electric  field  inside  and  outside. the  slab.  Substitu¬ 
ting  these  components  into  the  boundary  conditions  then  leads  to  a 
dispersion  relation,  whose  solutions  are  given  for  the  cases  shown  in 
fig.  2 • 
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The  following  symbols  are  used. 


E  (*f)  •  cross  sectional  distribution  of  the  longitudinal  electric 
!*  field  inside  the  plasma. 

(u. ^  i  cross  sectional  distribution  of  the  longitudinal  electric 
field  outside  the  plasma. 

E  E  s  amplitude  constants  of  Ei  (u.). 

o|(  Zf> 

i  longitudinal  propagation  constant 
LL  zz.  I" b  X 
2.  cL  :  slab  thickness 
2.  D  :  distance  between  plates. 

s  value  of  x  at  which  the  electron  density  is  zero . 

&  I'b  cl 

A  r  ftD 

X  =  pi 

fM  :  cross  sectional  distribution  of  the  electron  density. 

(_  :  second  order  differential  linear  operator. 

jJi,  s  parameter  characterizing  the  cross-sectional  inhomogeneity, 
s  cyclotron  angular  frequency. 


:  plasma  angular  frequency  at  the  center  of  the  slab. 

k  = 

<*Z-  OlX/lJ 
x  =  i  _(t  -  c syio 

^0-1 -«/k 
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2.  'Derivation  of  the  equation  for  Ez  (lA,) . 

Starting  from  Maxwell's  equations  a  wave  equation  for 
E  =  (Ex  Ey  Ez)  is  obtained,  s 


Vt_VVE  +  o 


where  fc  is  the  dyadic  absolute  permittivity. 

If  quasi-neutrality  is  assumed, (dielectric  approximation) 


V-(t.E)  = 


and,  hence,  V'E  is  derived,  the  z  component  of  eq.(l),in  the  quasi¬ 
static  approximation,  can  be  written  in  the  form  : 

!£.  +  x  1L  Uk  _  Ax  E* ,  o  (3 

<D  x*  £,  3x  2  x  /  <5| 

where  <S  and  S  are  the  diagonal  components  of  £  and  an 

I  9 

exp  j  Z)  dependence  is  assumed. 

In  eq.(3)  the  second  term  is  due  to  the  inhomogeneity  in  the  cross- 
section  and  the  absence  of  any  component  of  the  magnetic  field  is 
due  to  the  quasi-static  approximation. 

If  the  normalized  variable  LL«(bxis  used  eq.(3)  can  be 


^  Ez  ^  d  ti  1)  jT 

OK,1  e,  D  u.  9  u.  <5,  z 
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The  following  distributions  in  the  cross  section  (plasma 


region)  are  now  assumed  s 

AH- 

^  bO*) 

■i  _  «z 

£>(a)= 

i  -  b  ( u.) 

a.  <0$. 

with  o(  =.  —  and 

bW=bo  j(u)=(^/^){(u) 

specifies  ±  particular  kind  of  density  distribution. 

If 

*1 

I 

II 

where 

IV  /  x 

e,M  *ni  4»M 

li 

r^o 

N> 

|  can  be  rewritten 

£,W  = 

i  b„  (£-  el) 

*('-«*) 

£,M  = 

l 

0 

_£) 

1 

Substituting  these  values  into  eq.(4)  the  following 
equation  is  obtained  : 

L(Ez)  =  (a*-ir)  |£»  +  a-  |E.  ^)(<-«*)Ez  =0 


where 


Tr  I- 


\  - 


T=  (  -  _L 

b. 

and  the .first  member  of  the  equation  has  been  designated  with  the 
operational  symbol  L  (Ez). 


(5) 
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3.  Solution  of  the  equation  f  or  ^(U.)  inside  the  plasma. 

A.  solution  of  eq.  (5)  in  the  plasma  region  can  he  obtained 
with  the  aid  of  the  Ritz  method  by  assuming  a  test  function  as  an 
admissible  analytical  expression  for  Ez(u). 

If  a  linear  combination  of  suitable  functions  is  chosen  for 


a  first  order  approximation  is  given  by 


(6). 


where 


S  - 1^>  d 


Introducing  the  expression  (6)  of  EZp(*A.)  in  eq.  (5)  the 
problem  is  reduced  to  a  "minimum"  problem  with  the  following  boundary 
conditions  : 


The  problem  can  be  mathematically  formulated  as  follows! 


and  if  the  degree  of  homogeneity  in  the  cross-  section  is  specified 
by  a  characteristic  parameter  W.  defined  as  follows  : 


foglio 
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eq.(8)  can  be  reduced  to  s 


4-  Solution  of  equation  for  Ez(  \x)  outside  the  plasma. 


Outside  the  plasma  eq.(5)  reduces  to  the  following  ordinary 
second  order  differential  equation  : 


<)  Ei  _  E*  =  o 
d  £ 


(10) 


By  imposing  the  boundary  condition  at  U.  s  /l  ,  the 
solution  of  eq.  (10)  is  : 


H  »  A  ( 


e  _  e  I 


(ii) 


where  A  is  an  arbitrary  amplitude  constant. 
5 .  Derivation  of  the  dispersion  relation. _ 


The  eigen  value  h  can  be  calculated  by  imposing  the 
boundary  conditions  (7)>  at 


When  this  is  done,  the  following  two  equations  are  obtained  i 


E  =  A 


(e  -e  ) 


(12) 


E  =  A_  (€ 


zSi,® 


-A 
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Substituting  these  expressions  into  eq..(9)>  * 


This  is  the  dispersion  relation  in  the  general  case  of 
>^X^©oand  S^-A^o©  .  Some  particular  cases  of  special  interest  can 
he  derived  from  this • they  are  shown  in  the  following  table. 


PHYSICAL 

STRUCTURE 


DENSITY 

DISTRIBUTION 


CUPl/E 

RESIGNATION 


a)  A  =  & 

( COMPLETELY 
FILLED  WAVEGUIDE) 


B)  A  -*  oo 

('SLAB-  IN 
FREE  SPACE ) 


=  1  ( parabolic) 


JX=oo  ( UNIFORM ) 


(U.=  1  ( parabolic) 


1^/jL  Aoo  (PARABOLIC) 
(FOR  EX. yU*=  z) 


mm 


JA=°o  (  UNIFORM ) 


Particular  cases. 

In  the  a)  case  if  we  insert  Z\  -  c>  in  eq.  13)  the  R.H.S 
vanishes  and  the  following  equation  is  left  : 


or  the  alternative  form  : 


If  yU.«|  the  above  equation  reduces  to 

hl:_YL  l.-HQ-*) 

<  -  b0- 


(15) 


(16) 


which  is  the  dispersion  equation  of  a  completely  filled  waveguide  with 
a  parabolic  electron  density  distribution  in  the  cross-section. 


A  plot  of  eq.  (l6)  in  the  plane  v$  (&|  yields  the 

Wf  bo 

curves  designated  by  3^  .  (Each  curve  has  been  calculated  for  a  fixed 

a  z. 

value  of  the  parameter  O.;  _Wfc> ( see  tables  I  and  I'). 

Wt  K 


If  yit  =  6<3eq.  (15)  becomes 


If  ot-O  (isotropic  plasma) 

b\-l 

£ 

there  is  no  propagation. 


r  O 


(17) 
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Ifa(^o(gyrotropic  plasma) 

l  -  ^  b„  -  / 


(18) 


which  is  the  dispersion  relation  of  a  completely  filled  waveguide  with 
a  uniform  distribution  of  the  electron  density  in  the  cross-section. 

If  eq.  (l8)  is  plotted  in  the  U)  _(c>j  plane, the  curves  3_  are 

I  OO 

obtained.  (Tables  I,  I'). 

In  the  b)  case  if  we  let  o°  in  eq.  (13),  this  eq. 

reduces  to+  : 


--Pii  '-Kv-fc 


(19) 


ft  - 


-2. 


which  in  terms  of  the  normalized  to  ,  has  the  form; 


When  tends  to  infinity  our  physical  structux'e  reduces  to  a  slab 
of  plasma  placed  in  an  unlimited  space  .  In  this  case  th  (A-S)  *=  / 
and  the  absolute  value  of  S  must  be  considered,  if  reciprocity  is 
to  be  respected. 

Prom  a  mathematical  stand  point  it  is  to  be  noticed  that  when  ^> <  0 
&v*v 


the  limit 


tends  to 


!ria 


D- 

in  the  same  manner  as  when 


|5>0. 
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For  OL=i- 


(zj£  - ) 

-%e”  w  ] 

*3 

1 

-1 

i 

H4| 

L _ 

\  '6  lo^j 

j 

[3  »5  ^  j 

C20) 


■ - 1 

l 

X 

l 

\ 

jSj  +  2.o/ 

i^3  -  !/3 

L  /p-i)j 

1 

_  2^-f 

a  O 


For  Q_  =  2. 


>  .-2. 


2.  M-  _  JL  _  2.  m?  uj' 

16  lo«,  ~rit 


’  + 


2.  2  ^ 
yi^ _ L  M-  to2. 


(20') 


2 

}M  -  f 


+ 


-  2 

u) 


/Vb~  V6  _  / 
5 


uf-2. 


=  o 


If  yM,  »  1  eq.  20,  20'  reduce  to  : 


/  JiL  _  UJZ]c>  4.  [jL  _  j_  to^jSL  uj2 

(lob  15  J  (lb  3>  /'  I  5  -.2- 


14.  _i\=o 

^w2--  (  3 


(21) 


(21') 


Eqs.  (21  )  are  the  dispersion  relations  for  a  slab  of  plasma  in  free 
space  with  a  parabolic  electron  density  distribution  in  the  cross 
section,  (curves  b^  tables  II,  II'). 

If  U/  «  oo  ,  eqg  20  reduce  to  : 


_2_ 

5 


,  zl 
■(_ul 


-Z 

to 


b%  (i-w*)  w^i[J,|  +  D2  5 -it2  „ 


=  O 


uT-2 


(22) 


(22') 
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which  are  the  dispersion  relations  for  a  slab  of  plasma  in  free 
space  with  a  uniform  electron  density  distribution  in  the  cross- 
section.  (Curves  b  ;  tables  II,  II').  Dispersion  curves,  for  the  case 

i  .  . 

P-  2  have  also  been  calculated  and  plotted  (curves  b2). 

7 .  Comparison  between  the  approximate  solution  and  the  exact  solution. 

A  direct  analysis  of  the  admissible  propagation  modes  in 
a  slab  of  uniform  plasma  in  free  space  shows  that  a  cosinus  distribu¬ 
tion  of  Ez  in  the  cross-section  is  possible  when  a  quasi  static  E 
mode  is  considered. 

Imposition  of  boundary  conditions  yields  a  characteristic 
equation  and  a  dispersion  relation  which  is 


where  the  modulus  of  [b|  has  been  chosen  in  accordance  with  recipro¬ 
city  principle. 

_  2. 

In  terms  of  the  new  normalized  variable  UJ  eq.  (23)  can 
be  rewritten  s 


This  equation  can  be  solved  numerically  and  dispersion  curves 


can  be  derived  and  compared  with  curve  b^  .(Table  IV,  IV'). (  A  subscript 
"e"  is  used  to  indicate  "exact  solution"  curves). 
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CONCLUSIONS 


Application  of  a  variational  method  to  a  quasi-static  problem, 
has  been  made.  The  choice  of  the  test  function,  based  on  the  knowledge 
of  the  field  distribution  in  the  uniform  plasma  case,  represents  a  first 
order  mathematical  approximation  capable  of  describing  a  mode  with  the 
lowest  order  transversal  distribution. 


An  inspection  of  the  dispersion  curves,  shows  that,  within  the 
above  limits  : 

1  .  -  passing  from  cases  3-^and  bj(uniform  electron  density  distribution) 

to  cases  3l  and  (parabolic  electron  density  distribution)  the  upper 
band  shrinks. 


-  For  the  cL(  curve,  when  the  parameter  cl  equals  the  critical 

value  0,456  the  upper  band  shrinks  to  zero  and  if  0,4.0.  ,propai 

gation  occurs  only  in  the  pass  band  Cl 

For  the  b  curves  (slab  of  plasma  in  f.ee  space)  critical  values  of 

z 

<X  have  been  computed  for  different  yd,  and  the  results  plotted  in 
fig.  3. 
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-  The  asymptotic  io  value,  equs.l  to  1  in  the  uniform  case,  in  the 
passage  from  3-^to  <3U(case,  decreases  whereas  the  asymptotic 
value  equal  to  cx  remains  unchanged.  This  is  expected,  since 
the  average  electron  density,  in  .  the  passage  from  oL^to  <3-(  case, 
decreases  as  well. 

-  Letting  the  metal  walls  move  away  from  the  plasma  boundary  does 

not  affect  the  behaviour  of  the  dispersion  curves  as  much  as  varying 
the  electron  distribution  in  the  cross-section . (Table  III  and  III'). 

-  A  comparative  exam  of  the  dispersion  curves  obtained  by  variational 
methods  and  by  exact  mathematical  analysis,  for  a  slab  of  plasma 

in  free  space,  shows  a  fairly  good  agreement  in  the  upper  band,  and 
an  excellent  agreement  in  the  lower  band  when  Uv>u).  (Table  IV'). 

When  IO, ^  u),  the  lower  branches  of  the  curves  show  a  different 

b 

behaviour.  (Table  IV). 

It  has,  in  fact,  been  demonstrated  that,  with  or  without  quasi  static 
approximation,  (see  ref.  5  )  the  lowest  order  mode  dispersion  curve 
approaches  the  asymptotic  value  J  +  whenever  whereas 

application  of  the  variational  method  provides  only  an  asymptotic 
value  for  = 
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